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We investigate soliton mobility in the disordered Ablowitz-Ladik (AL) model and the standard 
nonlinear Schrodinger (NLS) lattice with the help of an effective potential generalizing the Peierls- 
Nabarro potential. This potential results from deviation from integrability, resulting of randomness 
for the AL model, and of both randomness and lattice discreteness for the NLS lattice. Statistical 
properties of such a potential are analyzed, and it is shown how the soliton mobility is affected 
by its size. The usefulness of this effective potential in studying soliton dynamics is demonstrated 
numerically. Further we propose two ways the soliton transport in presence of disorder can be 
enhanced: one is to use specific realizations of randomness, and the other one is to consider a 
specific soliton pair. 
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Introduction —Discrete solitons (including 

breathers—in our discussion referred as the local¬ 
ized modes of soliton type but with temporarily 
oscillating profiles) for the nonlinear Schrodinger (NLS) 
type models on a lattice, have been theoretically and 
experimentally studied in various areas of ph^ics, 
e.g., in nonlinear optics and solid-state physics [ll-Q. 
These localized modes can exist in the interplay of the 
nonlinearity and discreteness, and their mobility is one 
of the most important aspects, for considering the mass 
and energy transport in a lattice. One effective path to 
physically understand mobility of the discrete soliton 
is using the description in terms of the Peierls-Nabarro 
(PN) potentials [HI, Q. It could be understood as a 
type of effective potential, resulting of deviation from 
integrability (we will obtain similar deviation results of 
randomness for the disordered AL model, and of both 
randomness and discreteness for the disordered NLS 
lattice), with its amplitude being viewed as a minimum 
barrier that the soliton should overcome to propagate in 
the lattice. By virtue of such a concept, a series of works 
investigated mobility of discrete solitons in various NLS 
lattices [6|-[8|. 

One practical method to calculate the effective poten¬ 
tial is employing a collective coordinate approximation 
[HI, 0 , with the assumption that the discrete soliton prop¬ 
agates very slowly and adiabatically. In this potential the 
soliton can be considered as an effective particle. Such 
adiabatic approach provides a good approximate picture 
when the radiation is negligible. 

In the present paper, we will study solitons for the AL 
model and the standard NLS lattice in presence of dis¬ 
order using similar methods. Our investigations are rel¬ 
evant for experiments using random media in optics [9| 
and optical speckle potentials in Bose-Einstein conden- 
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sates M- On the other hand, nonlinearity and disorder 
are found to play important roles in existence and sta¬ 
bility of the localized excitation in a trimer model ini- 
We will analyze statistical properties of the generalized 
PN effective potential, and show how these affect the 
soliton mobility. We will mainly focus on the relatively 
short-time-scale soliton behavior, with the lattice soliton 
initially injected with a small amount of kinetic energy 
(moving slowly), and emitting negligible radiation. Un¬ 
der such situations, certain basic features related to the 
PN effective potential are believed to be revealed. 

Soliton mobility in Ablowitz-Ladik lattice with 
disorder —We start from the integrable Ablowitz-Ladik 
(AL) model [T2|, with a random potential term, defined 
by 

~ {fPn — l “k V^n+l)!! “k iV^nl ) “k , (1) 

where is the wavefunction at site n and time t, while 
Sn is the normally distributed random potential (uncor¬ 
related) with zero mean value and standard variance a, 
that is, {sn) = 0 and {snSn') = (J^S{n — n'). Integrability 
of the AL model in absence of the potential leads its PN 
potential barrier to vanish [H| (as well as vanishing of the 
resonances that generate radiation 0), and admits an 
exact mobile soliton solution at arbitrary velocity, writ¬ 
ten as 

. . . sinh(/i) 

Dn{t) = - -- -- exp[ik[n - x) + ia\ , (2) 

cosh[/i(n — x)\ 

where (without disorder) the time-dependent parameters 
X and a can be expressed as T = 2 sm{k) and a = 

2[cosh(/i) cos{k) + ^ sinh(/i) sin(/c)]. We now assume the 
random potential is weak and the soliton velocity is slow, 
then apply the adiabatic approximation MM to derive 
the evolution equations for the soliton parameters x and 
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k with disorder (assuming /i = 0) 
2 sinh(/i) 


X = 


+ 00 

^ = - E 


■ sm{k) , 


<^(^) = - coshlMt+tVcTshK-i)] • The force /(C) also has 
zero mean value, and its correlation function r/ can be 
(3a) derived due to its linear property, 


Sn sinh^(/i) tanh[/i(n — x)] 
cosh[//(n 1 — x)] cosh[/i(n — 1 — x)] 

(3b) 


+ 00 


r/(c,e') = (/(e)/(r)) = ^" E Hn-o^n-a- (s) 


Using the Poisson summation formula and residue theo- 
For the system ([3j), the effective potential, regarded as rem, we can evaluate the sums in m, and obtain 
the PN potential, can be written as 


u{x )=- r /(c)dc, 

Jo 


( 4 ) 




+ 00 


lHo(M,AC) + 2^fR,(M,C,AC) 


S = 1 


, (6) 


with the equivalent force /(^) = '^en4>{Ti — f), where where (A^ = f — ^'), 

I 

fHo =sinh^(/i){2A^coth(/iA^)csch[/i(A^ — l)]csch[/i(A^ + 1)] — (A^ — 2) coth[/i(A^ — l)]csch[/i(A^ — 2)]csch(/iA^) 


-(A^ + 2) coth[/i(A^ + l)]csch[/i(A^ + 2)]csch(/iA^)} , 


(7a) 


93s ^^^ csch cos[7r5(2^ — A^)]{2 coth(/iA^)csch[/i(A^ + l)]csch[/i(A^ — 1)] sin(7r5A^) 

li \ ji J 

— coth[/i(A^ — l)]csch(/iA^)csch[/i(A^ — 2)] sin[7r5(A^ — 2)] 

— coth[/i(A^ + l)]csch(/iA^)csch[/i(A^ + 2)] sin[7rs(A^ + 2)]} . (7b) 


Eqs. ([6]) and ([7|) apparently show that, due to discrete- not only on A^) [T5|. With (|6]), we derive the variance 

ness, /(^) is a nonstationary random process (depending of /(^) as 

I 


q f\ ^ +CO / 2 \ 

^ 0= 3 coth(/i) — tanh(/i)-^ ^ 5 csch f j cos(27rs^) . (8) 


Generally speaking, when // <C 1, f{f) can be ap¬ 
proximately seen as a stationary random process, with 
cTj/cr^ ^ 3coth(/i) — tanh(/i) — ^ (when // < 0.3, it 
accords well with the continuous limit result = 

~ when 1 < M ^ 3, /(C) is a non- 

stationary random process with a periodic variance (the 
terms for s ^ 2 are neglected as small terms) ~ 

3coth(/i) — tanh(/i) — ^ — ^-csch cos(27r^); when 

M ^ 3, we may have to consider overlap of the terms 
with 5^2. We numerically calculate the force /(^) for 
a large number of realizations of the random potential 
to derive its variance, and make a comparison with (©. 
Good agreement can be seen in Fig.^a), for two different 
regimes. 

With the statistical property of the force /(O? we next 
consider the effective potential U{x). Apparently, it is a 
nonstationary random process with zero mean value, and 


its correlation function can be derived as 

ru{x,x') = {U{x)U(x')) = j [ rf{^,C)d4dC ■ (9) 

Jo Jo 

We can numerically integrate (|9]) with appropriate trun¬ 
cation of ([6]) according to the value of /i, and further 
obtain the variance of U(x) as 

g}j{x)Ig^ = ru(x' -)> x)/(T^ . (10) 

Three typical examples of m, compared to numerical 
results found sampling various realizations using (j4]), are 
shown in Fig. [Hb). We can see that, after an increase in 
short distance, au/o- approaches a periodic type function 
(for small enough //, the periodicity can be neglected, 
e.g., /i = 0.5) with a nonzero mean value (averaged in 
finite periods along x). We denote such mean value as 
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I where is theoretically defined as 

= lim ^ [ cFu{x)dx , (11) 

L^oo L Jq 

and further compute that the amplitude of the periodic 
function is not more than 6% of even up to /i = 6. 

Thus, can be employed, in the sense of statistics, as 
the PN-type potential barrier that relates the soliton mo¬ 
bility in presence of disorder. cr[7^^/cr as a function of /i, 
is presented in Fig. [U^c), as well as the comparison with 
numerical results for multiple realizations of the random 
potentials. Generally speaking, the smaller soliton has 
a larger mobility due to its smaller effective potential 
barrier, however, there exist approximately two regimes: 
one is the sharp slope regime (/i < 1) where fa de¬ 
creases quickly with /i decreasing; the other one is the 
flat slope regime (/i > 3) where a^^ fa decreases consid¬ 
erably slowly with /i decreasing. 

To check such mobility, we directly proposed numeri¬ 
cal simulation of Eq. m, with the soliton solution ([2]) as 
the initial condition. The soliton behavior was observed 
for more than 100 realizations of random potentials with 
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FIG. 1. (Color online) Statistical properties of the equiva¬ 
lent force /(^) and effective potential U{x). The blue solid 
lines represent the results from jS]) or CQl), and the red circles 
represent the statistically simulated results. In simulation, 
we discretize the coordinate with = 0.05, and the vari¬ 
ance is computed for 10^ realizations of random potentials, 
(a) Comparison between the statistically simulated variance 
cTj/cr^ and (|8|) for /x = 0.5, 2.0, and 4.0 {a = 0.02). (b) Com¬ 
parison between the statistically simulated variance afjfa^ 
and (flOl) . (c) fa of (fTTl) as a function of ji. 


a = 0.01, and the soliton was given a small initial ve¬ 
locity Vo = 0.3 for each realization. Two typical results, 
with the soliton center of mass Xc as function of t, are 
presented in Figs. EJa) and (b). We can see that, no 
matter whether the soliton decelerates or accelerates (on 
average) for some time, the larger soliton is reflected or 
trapped by the randomness stronger than the small soli¬ 
ton. In principle, the soliton reflects near some large 
effective potential barrier, or gets trapped between two 
large barriers. Before the first reflection, whether deceler¬ 
ation or acceleration occurs, is depending on the details of 
the random potential. Generally, we can approximately 
use a time difference to identify such deceleration (accel¬ 
eration) stage, defined as 


fX 

t{x) = dx 

Jo 


2 sinh(/i) sm{k) 


1 

'^0 


( 12 ) 


where cos{k) = 2 sirih(/x) ^ cos(/co), as derived from 
the conservation of the effective total energy H = 
— 2sin^(^) cos{k) + U{x). Thus, if r < 0, there exists an 
averaged acceleration process within [0,x] (before first 
reflection); while r > 0, it is an averaged deceleration 
process. One interesting idea is to use the region where 
randomness accelerates the soliton in finite time. For in¬ 
stance, as seen in Fig.[2fb), for relatively large soliton, we 
may choose the random potential sections before the first 
reflection point, and arrange them periodically, to realize 
the mass transport in one direction for some time. 

Another aspect is to consider a pair of solitons, where 
one is large and the other one is small, with the same 
initial velocity and a short separation distance between 
them (may be partially overlapped). In order to enhance 
transportation of the large soliton, we arrange the small 
one to follow the large one, and give them an initial phase 
difference tt that provides a repulsive interaction between 
both solitons. When the large soliton encounters a large 
effective potential barrier that would obstruct its mo¬ 
bility, while the small one might be still well propagat¬ 
ing due to its larger mobility, and is expected to push 
the large one to help it overcome the barrier. A simi¬ 
lar mechanism to trigger a migration by a low-amplitude 



FIG. 2. (Color online) Soliton center of mass Xc as a function 
of time t for two typical realizations of random potentials 
(cr = 0.01 and vq = 0.3), found from ([Tj). 
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solitary wave that collides with an excitation has been 
described for a pure NLS lattice [T6|. As an example, 
such idea is realized and illustrated in Fig. [3] for a spe¬ 
cific realization of randomness. Fig.[3fa) presents a large 
single soliton (/i = 2) trapped by the random potential, 
while in Fig. [3l^b), we add the auxiliary small soliton 
(/i = 0.5) with the separation distance of 5. It shows 
that the large soliton can now propagate without reflec¬ 
tion or trapping for some time. The interaction of the 
small soliton (/i = 0.5), with the large one results in the 
motion of the large soliton in the same direction with¬ 
out reflection. Similar results were found also for other 
realizations of the disorder. 

For the large soliton, the randomness-induced radia¬ 
tion is still very small for considerably long time, and 
the soliton behavior can be well described by the effective 
potential [ a typical comparison is shown in Fig. ija) up 
to i = 10^ ]. Generally speaking, the large soliton is apt 
to be trapped by the randomness after some sequences 
of deceleration and acceleration periods (the soliton be¬ 
havior was observed in the long-time simulation for more 
than 10 realizations of the random potential). For the 
small soliton, if the randomness is weak enough, the ef¬ 
fective potential also gives a good approximation for long 
time. In Fig. IH^b), we decrease the strength of the ran¬ 
dom potential from the left to the right panels, and show 
that the effective potential approach agrees better with 
the numerical solutions of Eq. m for weaker random¬ 
ness. On the other hand, if the randomness is relatively 
strong, the small soliton was observed to continue radi¬ 
ating its mass and kinetic energy during the long-time 
propagation, and even exhibit visible deformation on its 
profile. Such a case may be analyzed using the method of 
modulation equation [^, which is not within the scope 
of this paper. 

Behavior for the standard NLS lattice with dis¬ 
order —We study the standard NLS lattice 

ii’n = -{'tpn-l + V’n+l) “ i^lV’nPV’n + ^nV’n • (13) 


We consider this equation as the AL model with a per¬ 
turbation term on the RHS, Rn = |'0nP('0n+i + V^n-i ~ 
z/^/in) and use the adiabatic approximation to de¬ 

rive the similar generalized PN effective potential when 
fi is not too large (/r ^ 3), 

U{x) =- j ^ cos(27rx) + t/(x) , (14) 

sinh f ^ j 

where /i = 0 and x 2 sinh(/i) tanh(/i) sin(/c)//i^, and 
U{x) is the same randomness generated effective poten¬ 
tial 0. In the limit /i ^ 0, the soliton is wide, therefore 
the discreteness of the lattice is not important, and the 
result of the continuous integrable model is approached. 
Naturally, we may use the following parameter to approx¬ 
imately decide which factor dominates, 

(^) 27r^z/sinh^(/i) a ^ 

^ 27r2sinh2(^) 

(15) 

where A(/i) denotes the curve presented in Fig. [TJc). For 
the typical parameters cr = 0.01 and = 1, when /i ^ 1, 
^ 1, the randomness dominates the soliton behavior; 
when /i 1, K ~ 2.6, these two potentials are of the 
same magnitude; when ja ^ k the influence of 

randomness can be ignored. And, if the soliton is not too 
large, influence of the randomness on the soliton mobility, 
with regard to its size, is similar to that found for the 
model ([I]). 

Here we make some comments: in fact, if the soliton 
is very small, it is easy to emit relatively strong radia¬ 
tion induced by the randomness, even to be greatly de¬ 
stroyed on its profile. Such a condition may be remark¬ 
ably out of the adiabatic approximation. If the soliton 
is very large, its mobility becomes much smaller, due to 
the large potential barrier of discreteness. On the other 
hand, strictly speaking, the soliton discussed here should 



(a) (b) 

FIG. 3. (Color online) Intensity map of the soliton propa¬ 
gation in a realization of random potential with a — 0.01. 
(a) A single soliton with /x = 2 and vo = 0.3. (b) A pair of 
solitons with fii = 2, /X 2 = 0.5, and vq = 0.3. The initial 
separation distance Ax and phase difference Aa are 5 and tt, 
respectively. 





FIG. 4. (Color online) Comparison of the soliton trajec¬ 
tory Xc obtained by numerically integrating Eq. ([1]) (red solid 
line) and effective potential approach using U{x) of (|3|) (green 
dashed line), (a) The random potential e used is of cr = 0.01, 
and the soliton parameters are /x = 1.0 and vq = 0.3. (b) 
Left Panel: the random potential is s' = 0.5^; Right Panel: 
s' = 0.26. The soliton parameters are /x = 0.5 and vq = 0.3. 
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be replaced by the breather (/i 7 ^ 0 ), however, for some 
parameters, our study may generally and effectively pro¬ 
vides a physical description of the lattice soliton mobility 
with disorder. As an example, we give a set of simu¬ 
lation results in Fig. [5l We know that, Eq. (HSI with 
5 n = 0 has no exact mobile soliton (breather) solutions, 
however, we could start from a sech type soliton, after 
initially emitting a small part of radiation, to numeri¬ 
cally generate an approximate breather solution that can 
propagate for considerably long time [l 8 |. InFig. [5fa), 
we generate such two breathers, one is small (the upper 
panel), the other one is large (the lower panel), with al¬ 
most the same small velocity. Then, we add a realization 
of the random potential, and show the simulation result 
in Fig. mb). Apparently, we can see that, the mobility 
of the large breather is obstructed by the randomness, 
while the small one is nearly not affected, which is the 
similar feature as shown in Fig. [2l 

Conclusions —With help of the generalized PN 
effective potential, we investigated mobility of lattice 
solitons in presence of disorder. We analyzed methods 
to enhance the mobility of solitons of the AL model ([T]) 
in a weak random potential. In some situations (weak 
randomness and short time), these can be considered as 
particles moving in the effective potential U{x) of (|1]). 
We find two ways to enhance mobility: (a) introducing 
a random potential that is engineered so that intervals 
of the random potential where acceleration takes place 
appear one after the other; (b) a small soliton is arranged 
to push a large one. The effective potential acting on 
the soliton results of the deviation from integrability. 
For the AL model o that is integrable in absence of 
randomness, the potential results of the randomness. 
For the standard NFS lattice ([T3|) . the deviation from in¬ 
tegrability results of both randomness and discreteness. 


as it is integrable in the continuum limit in absence of 
randomness. The effect of the deviations results in the 



(a) (b) 

FIG. 5. (Color online) (a) Intensity \'ipn\‘^ for two generated 
breathers of the NLS lattice ([13)) {sn = 0) with almost the 
same velocity v ~ 0.30. The initial conditions are chosen as 
([ 2 ]) with /i = 0.60 and k = 0.150 for the upper panel, while 
fi — 0.95 and k = 0.178 for the lower panel {n = 0.9). (b) 
The two breathers are propagating in presence of the random 
potential with a = 0.018. 

potential ([T4)) that plays the role of U{x) of (j4]) found 
for the AL model o. Therefore similar dynamics of 
solitons is expected in the regime where the effective 
potential description is valid. 


Z.-Y. S. acknowledges the support in part at the Tech- 
nion by a fellowship of the Israel Council for Higher Ed¬ 
ucation. This work was partly supported by the Israel 
Science Foundation (ISF-1028), by the US-Israel Bina¬ 
tional Science Foundation (BSF-2010132), by the USA 
National Science Foundation (NSF DMS 1201394) and 
by the Shlomo Kaplansky academic chair. 


[1] F. Lederer, G. 1. Stegeman, D. N. Christodoulides, G. 
Assanto, M. Segev, and Y. Silberberg, Phys. Rep. 463, 
1 (2008). 

[2] S. Flach and A. V. Gorbach, Phys. Rep. 467 , 1 (2008). 

[3] S. Aubry, Physica D 216 , 1 (2006). 

[4] P. G. Kevrekidis, The Discrete Nonlinear Schrodinger 
Equation^ Springer-Verlag (Berlin Heidelberg) (2009). 

[5] Y. S. Kivshar and D. K. Campbell, Phys. Rev. E 48, 
3077 (1993). 

[6] M. Johansson and P. Jason, arXiv: 1412.0994 (2014). 

[7] L. Hadzievski, A. Maluckov, M. Stepic, and D. Kip, Phys. 
Rev. Lett. 93 , 033901 (2004). 

[8] U. Naether, R. A. Vicencio, and M. Stepic, Opt. Lett. 
36 , 1467 (2011). 

[9] T. Schwartz, G. Bartal, S. Fishman, and M. Segev, Na¬ 
ture (London) 446, 52 (2007); Y. Lahini, A. Avidan, F. 
Pozzi, M. Sorel, R. Morandotti, D. N. Christodoulides, 
and Y. Silberberg, Phys. Rev. Lett. 100 , 013906 (2008). 

[10] J. E. Lye, L. Fallani, M. Modugno, D. S. Wiersma, 
C. Fort, and M. Inguscio, Phys. Rev. Lett. 95 , 070401 
(2005); D. Clement, A. F. Varon, M. Hugbart, J. A. Ret- 
ter, P. Bouyer, L. Sanchez-Palencia, D. M. Gangardt, 


G. V. Shlyapnikov, and A. Aspect, Phys. Rev. Lett. 95 , 
170409 (2005); C. Fort, L. Fallani, V. Guarrera, J. E. Lye, 
M. Modugno, D. S. Wiersma, and M. Inguscio, Phys. 
Rev. Lett. 95 , 170410 (2005). 

[11] X.-D. Bai and J.-K. Xue, Phys. Rev. E 86, 066605 (2012). 

[12] M. J. Ablowitz and J. F. Ladik, J. Math. Phys. 16 , 598 
(1975); 17 , 1011 (1976). 

[13] A. A. Vakhnenko and Yu. B. Gaididei, Tero. Mat. Fiz. 
68, 350 (1986) [Theor. Math. Phys. 68, 873 (1987)]. 

[14] D. Cai, A. R. Bishop, and N. Gr0nbech-Jensen, Phys. 
Rev. E 53, 4131 (1996). 

[15] A. Papoulis and S. U. Pillai, Probability, Random Vari¬ 
ables and Stochastic Processes, 4th edition, McGraw-Hill 
Europe (2002). 

[16] B. Rumpf, Phys. Rev. E 70 , 016609 (2004). 

[17] A. Soffer, Proc. Int. Congress of Mathematicians 
(Madrid, Spain, 2006), 459 (2006); A. Soffer and X. Zhao, 
J. Phys. A 48, 135201 (2015). 

[18] R. Franzosi, R. Livi, G.L. Oppo, and A. Politi, Nonlin¬ 
earity 24 , R89 (2011). 





